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Abstract 


1 Introduction 


We present a detailed analysis of the tempera¬ 
ture statistics in an oceanographic observational 
dataset. The data are collected using a moored 
array of thermistors, 100 m tall and starting 5 m 
above the bottom, deployed during four months 
above the slopes of a Seamount in the north¬ 
eastern Atlantic Ocean. Turbulence at this loca¬ 
tion is strongly affected by the semidiurnal tidal 
wave. Mean stratification is stable in the entire 
dataset. We compute structure functions, of or¬ 
der up to 10, of the distributions of temperature 
increments. Strong intermittency is observed, in 
particular, during the downslope phase of the tide, 
and farther from the solid bottom. In the lower 
half of the mooring during the upslope phase, the 
temperature statistics are consistent with those of 
a passive scalar. In the upper half of the moor¬ 
ing, the temperature statistics deviate from those 
of a passive scalar, and evidence of turbulent con¬ 
vective activity is found. The downslope phase is 
generally thought to be more shear-dominated, but 
our results suggest on the other hand that convec¬ 
tive activity is present. High-order moments also 
show that the turbulence scaling behaviour breaks 
at a well-defined scale (of the order of the buoy¬ 
ancy length scale), which is however dependent on 
the flow state (tidal phase, height above the bot¬ 
tom). At larger scales, wave motions are dominant. 
We suggest that our results could provide an impor¬ 
tant reference for laboratory and numerical studies 
of mixing in geophysical flows. 
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Turbulence above sloping ocean boundaries has 
been the subject of several studies since the work 


of Munk (1966), who hypothesised that enhanced 


turbulence close to the solid boundary plays an im¬ 
portant role in maintaining the ocean stratified in 
density from the surface to the bottom. The dy¬ 
namics of the benthic boundary layer above a slope 


has been the subject of several studies (e.g. Garrett 



Thorpe et al.l 

1990| Garrettl 19911 |van Haren 

et al. 

|1994 Alford 

and Pinkel[ 2000 |Slinn and 


Levine[ 2003 Gayen and Sarka^ 20lia|b ). Tur¬ 


bulence close to sloping boundaries is enhanced by 
the breaking of internal waves (see the review of 
Lamb 2014) which, from a broader point of view. 


dissipates tidal energy (Munk and Wunsch 1998) 
and, to a lesser extent, eddy kinetic energy (Zhai 


et ah 


2010 ) 


Here, we consider turbulence above a deep-ocean 
sloping boundary from the particular perspective 
of scalar (temperature) statistics. The interest in 
the scalar statistics embedded in a turbulent flow 
is manifold. The scalar field is highly intermittent, 
and its statistical moments bear the marks of this 
intermittency. In relation to this, the study of the 
statistics of a scalar is the basis for understanding 
the dynamics driving its mixing, and for reliably 
estimating and modelling mixing itself. The statis¬ 
tics of scalar tracers in a turbulent flow is a widely 
studied topic in the laboratory as well as numeri¬ 


cally and analytically (see the reviews of Shraiman 
and Siggia 2000 War haft [ |2000 ); most often the 
focus is on grid-generated turbulence. 

On the other hand, intermittency of scalar fields 
in strongly stratified turbulent flows, where the 
buoyancy force is expected to be dominant, has re- 
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ceived limited attention (e.g. |BrethQuwer and LincT 


borg 2008 using numerical simulations). While 


we are unaware of laboratory studies on this topic, 


Seuront et ah (1999) discussed temperature fluc¬ 


tuations in the surface ocean boundary layer in 
terms of multifractal intermittency. [Costa Frola 


probability density function (pdf) of temperature 
increments in the atmospheric boundary layer, dis¬ 
cussing it in terms of a bifractal intermittency 
model. 

We must stress that observations in the ocean 
and atmosphere often include some temperature 
statistics, but these are virtually always limited to 
second order, i.e. to spectra (see, for instance, [Riley 


sipation rate by using the Osborn and Cox (1972) 
model, and only second-order statistics are pre¬ 
sented. An exception is Thorpe et al. (1991), where 


those of the inertial range (Rorai et al. 2014). This 


tified. In the present observations, on the other 
hand, there is no obvious separation between break¬ 
ing waves (the main forcing) and turbulence. In¬ 
stabilities such as those produced in the laboratory 
are seldom recognised in the ocean, even if new 


et al. (2014) recently reported observations of the 


observations may change this (van Haren and Gos- 
tianx[ 2010 van Har^ [2013| [Smyth and Moum[ 


2012 and references therein). For this reason, the 


relevance of the conclusions drawn from the study 
of Kelvin-Helmholtz instabilities for the real ocean 
is still not clear (Thorpe 2012). In this context. 


and Lindborg 2008 and references therein). The 


oceanographic literature usually considers observa¬ 
tions of temperature at the dissipation scale with 
the aim of estimating turbulent diffusivity and dis¬ 


a possible link between the skewness of temperature 
increments and convection in benthic boundary lay¬ 
ers was discussed. Possibly more often, ocean tem¬ 
perature (or density) observations are collected at 
lower resolution, comparable to the Ozmidov scale 
(e.g. Alford and PinkeH 2000). 

This work presents a detailed study of the 
temperature statistics above a deep-ocean sloping 
boundary, focusing on the intermittency properties 
of the scalar. The results are discussed compar¬ 
ing them to laboratory experiments with both pas¬ 
sive (shear turbulence) and active (convective tur¬ 
bulence) scalars. By necessity, our observations 
provide an average view of fluid regions having 
different Reynolds number, different mean shear 
and different mean stratification. While this is 
to some extent a limitation, it should be kept in 
mind that geophysical turbulence is intrinsically 
sporadic i.e. intermittent also at scales larger than 


fact must be taken into account when studying the 
scalar statistics in a stratified flow, and is naturally 
included in real-ocean observations. Furthermore, 
laboratory and numerical experiments mainly focus 
on grid-generated turbulence or on the evolution 
of a particular kind of instability, usually Kelvin- 
Helmholtz. A unique forcing scale (grid length, 
most unstable wavenumber) can usually be iden¬ 


the main aim of this work is to help bridge the 
gap between laboratory and geophysical flows, by 
analysing the observations in a way that is familiar 
to the fluid dynamics community (pdf’s, moments, 
structure functions), and thus could serve as a ref¬ 
erence. 

The data used in this work have been collected 
using a moored array of thermistors deployed above 
the slopes of Seamount Josephine, in the north¬ 
eastern Atlantic Ocean. Motions at this loca¬ 
tion are dominated by the tide. The thermistors 
used, built in-house at the Royal Netherlands In¬ 
stitute for Sea Research (NIOZ), measure tempera¬ 
ture variations at a relatively high temporal resolu¬ 
tion (1 s) for long periods of time (several months). 
While their temporal resolution is generally insuf¬ 
ficient to resolve the complete inertial range of tur¬ 
bulence (nor to resolve the dissipation scale), they 
do capture part of the turbulence cascade. In this 
work, we concentrate on time-averaged statistics, 
distinguishing between the upslope and the downs- 
lope phases of the tide, and distinguishing differ¬ 
ent parts of the moored array, i.e. different depth 
ranges. On the other hand, we do not consider time 
modulation of turbulence properties other than the 
tidal one, i.e. we do not consider the modulation 
by inertial waves, mesoscale eddies or other low- 
frequency processes. Time averages are a necessary 
starting point, enabling better-converged statistics, 
and thus enabling higher-order statistics to be reli¬ 
ably computed. 

We will consider almost exclusively fluctuations 
and increments of temperature within the horizon¬ 
tal plane. This is motivated by the fact that the 
resolution of our data is substantially better in 
the horizontal direction (0.2 m by using Taylor’s 
hypothesis; see section [2.2[ ) than in the vertical 
(0.7 m). The difference in resolution is particularly 
relevant if the anisotropy of turbulence in a strongly 
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stratified environment is taken into account, with 
length scales much smaller along the vertical di¬ 
rection (approximately parallel to the mean den¬ 
sity gradient) than in the horizontal (approximately 
perpendicular to the mean density gradient) |Billant| 
and Chomaz| (2001). 

The paper is organised as follows. The tempera¬ 
ture data used throughout the paper are described 
in section [2T| Velocity observations are briefly de¬ 
scribed in section [ 2 I 2 I The statistical observables 


considered are concisely presented in section 2.3 


The results are presented in section discussion 
and conclusions follow in section 01 


2 Methods 

2.1 Temperature data 

The data used in this work are obtained from a 
mooring deployed between spring and late summer 
2013 in the north-eastern Atlantic Ocean on the 
slopes of Seamount Josephine (see table for de¬ 
tails). The mooring has 144 “NIOZ4” thermistors 
(an evolved version of the ones described in 


Haren et al. 


2QQ9| taped on a nylon-coated steel 


cable at intervals of 0.7 m. The thermistors sam¬ 
pled temperature at a rate of 1 Hz with a preci¬ 
sion better than 1 mK. The time series from each 
thermistor is approximately 10^ s long, with the 
entire dataset containing approximately 10^ mea¬ 
surements. The cable was attached to a 500 kg 
ballast weight at the bottom, and to an ellipti¬ 
cal buoy at the top. The high net buoyancy (ap¬ 
proximately 300 kg) together with the small size of 
the sensors, limiting the drag on the cable, guaran¬ 
tee that measurements are very close to Eulerian, 


as checked also in similar previous moorings (van 


Haren and Gostiaux 2009 2012). The pressure and 


tilt sensors of the acoustic Doppler current profiler 
(ADCP) mounted in the top buoy, 210 m above the 
bottom, recorded movements smaller than 0.2 m in 
the vertical direction and 3 m in the horizontal one. 
All the thermistors performed satisfactorily for the 
whole deployment, with a noise level of approxi¬ 
mately 5 X 10“^ °C. Calibration is applied to the 
raw data and the drift in the response of the ther¬ 
mistor electronics, visible over periods longer than 
a few weeks, is compensated for. 

At the mooring location, temperature and ve¬ 


locity variance are dominated by the semidiurnal 
tidal period (in particular by the M 2 lunar period 
of 12.42 hour), as shown by the velocity spectrum 
in figure (left panel). Tidal variations are visible 
in the temperature record as successive cooling and 
warming trends. Throughout this paper, we con¬ 
sider various quantities computed separately for the 
cooling and warming tidal phases. For simplicity, 
we will usually refer to the cooling phase as “up- 
slope” tide and to the warming phase as “downs- 
lope” tide. The two phases are considered sepa¬ 
rately because statistical properties vary markedly 
between the two. The upslope and downslope 
phases are defined by the sign of the time deriva¬ 
tive (negative and positive, respectively), of the 
vertically averaged temperature signal, band pass- 
filtered at the semidiurnal frequency. This defini¬ 
tion produces two subsets of the data that have 
virtually the same size. Two snapshots are shown 
in figure The two panels show approximately 
30 min of data each, during downslope (labelled 
“down”) and upslope (labelled “up”) tidal phases, 
respectively. Figure also shows the four verti¬ 
cal segments of the mooring (A, B, C, D) consid¬ 
ered separately in the rest of the work, each ap¬ 
proximately 25 m high and including 36 thermis¬ 
tors (table 1^. From figure]^ we can see that, dur¬ 
ing both the upslope and the downslope phases, 
water is stratified at the mooring location. This 
is confirmed by table reporting mean, first and 
ninety-ninth percentiles of the buoyancy frequency 
(N) in the four segments of the mooring and dur¬ 
ing the two tidal phases separately. Stratification is 
computed from 25 s averages of temperature, with 
profiles vertically reordered to obtain a statically 
stable stratificatiorfl This observation is consis¬ 
tent with those by [van Haren and Gostiaux (2012), 
which extended down to within 0.5 m from the solid 
boundary, differently from those discussed here, 
having the deepest thermistor 5 m above the bot¬ 
tom. More details on the computation of buoy¬ 
ancy frequency using these data can be found in 
Gimatoribus et al. (2014). The mean buoyancy fre- 


^This procedure, used to obtain an estimate of the back¬ 
ground stratification, removes all instabilities from the verti¬ 
cal profiles. On the other hand, unstratified regions are not 
removed from profiles and would still produce vanishing val- 
ues of the buoy ancy frequency, if present. See [Winters and] 
jP’Asaro] ( |1996| > for a formal justification of the reordering 
procedure in the context of mixing estimation. 
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quency is slightly higher during the upslope tidal 
phase, with the exception of the lowest segment, 
segment A. 


2.2 Velocity data and Taylor’s hy¬ 
pothesis 

A Teledyne/RDI 75 kHz ADCP was mounted in 
the top buoy of the mooring, 210 m above the bot¬ 
tom. The ADCP was pointed downwards, resolving 
current and echo intensity observations every 5 m 
vertically, between 15 and 185 m above the bot¬ 
tom. Horizontal vector velocity is measured, but 
only the magnitude within the depth range of tem¬ 
perature measurements is considered here. We will 
refer to the magnitude of velocity simply as “veloc¬ 
ity” in the following, unless stated otherwise. The 
sampling interval was 600 s. 

Despite the limited resolution, velocity data can 
be used to transform the measured temperature 
from the temporal to the spatial domain, accord¬ 
ing to the Taylor’s hypothesis of frozen turbulence. 
Here, we must take into account that the mean ve¬ 
locity is not constant, following the approach of 
Pinton and Labbe (1994). The velocity held v is 


hrst decomposed into a mean component (u) and 
a huctuating one v' using a low-pass Butterworth 
hlter with stop band 1/cr = 3000 s. On average, 
this can be considered the forcing frequency, based 
on the frequency spectra of temperature huctua- 
tions (not shown). Approximately at this value, 
the spectral slope decreases moving towards higher 
frequencies, approaching the value of —5/3 typical 
of the inertial range of turbulence in both homo¬ 
geneous and stratihed environments dTennekes and 


Lumley[ [l972| Lindborg and Fedina |2009 ). A sim¬ 

ilar transition is observed in the velocity spectrum 
in hgure but in this case the turbulence range 
is mostly undetected and contaminated by noise 
(light blue line in the hgure). The pdf’s of u, {v) 
and v' are shown in hgure(right panel). Mean ve¬ 
locity is higher during the upslope phase of the tide, 
while the pdf’s of velocity huctuations are similar 
for the upslope and downslope phases. 

In order to transform the temperature time se¬ 
ries of each thermistor into the spatial domain, the 
mean velocity is averaged, along the vertical co¬ 
ordinate, within each of the four segments of the 
mooring. Since velocity is not constant in time, 
the spatial series thus obtained has a non-constant 


sampling rate. Linear interpolation is performed 
to a constant time step of 0.2 m, which reduces the 
size of the dataset by approximately a factor of 2. 

The results of this transformation have to be 
evaluated with care, since we cannot assess the 
goodness of the Taylor’s hypothesis, lacking in¬ 
formation on the turbulent velocity. This is par¬ 
ticularly true closer to the bottom, where turbu¬ 
lence and shear are stronger. Turbulence intensity 
((u'^) / (u)) computed from the ADCP data is 

much smaller than unity most of the time. While 
this gives support for the use of Taylor’s hypothe¬ 
sis, we do not emphasise this result since the ADCP 
resolution is insufficient (both in time and in the 
vertical) to detect most of the turbulence range, 
and noise contaminates the highest frequency range 
of ADCP measurement^ In the next subsections, 
the results are mostly based on the spatial (Taylor- 
transformed) data. The same quantities reported 
here have been computed also on the original (time 
series) data, with qualitatively similar results. Rel¬ 
atively short spatial separations are considered, 
which correspond in the time domain to time in¬ 
tervals no longer than approximately one hour. 


2.3 Statistical description 


A power spectrum can fully characterise a dataset 
only if the distribution of the data is approximately 
Gaussian. This is usually not the case for temper¬ 
ature (and scalars in general) in a turbulent flow, 
nor is it the case here. For this reason, our de¬ 
scription of the dataset is extended to higher-order 
moments. In doing so, we follow the example of 
the analyses of laboratory data available in the lit¬ 
erature, mostly concerning passive scalars in grid¬ 
generated turbulence; see for example the works by 


Thoroddsen and Van Atta (1992) and by Mydlarski 


and Warhaft (1998), and the reviews by Warhaft 


(2000) and Shraiman and Siggia (2000). We also 


make use of some of the concepts developed in |Zhou 
and Xia (2002) for convectively generated turbu¬ 


lence. We thus apply the now-classic methods of 
these works to the case of a strongly stratihed ocean 
environment, with the aim of contributing to the 
understanding of the latter. In fact, we hnd that 
the results from controlled environments are often 


^The presence of noise artificially increases the estimated 
turbulence intensity. 
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comparable to those from our observational data. 

We define here a common notation for the nor¬ 
malised moments of a distribution, defined for a 
generic variable x at order q by 


N (x) 


frequency (see table |^; qualitatively similar results 
are obtained using a higher stop-band frequency). 

If the quantity x in equation is an increment, 
whose value depends on the distance between two 
measurements (in either time or space), then fiq (x) 
depends on the same distance too, and is usually 
called a structure function. In particular, in section 
|3.3| we consider the temperature increments: 


with the average indicated by the angular brack¬ 
ets. With this notation, the third-order moment 
of X, that is the skewness, is written as /is (x) and 
the fourth-order moment, the flatness, as /i 4 (x), 
and both are defined by equation These two 
moments have received particular attention in the 
literature (War haft 2000), as a way to characterise 


ArO = 0{ro +r) - 6>(ro), 


( 2 ) 


intermittency and anisotropy (skewness). 

In section [3^ we consider the temperature fluc¬ 
tuations, 6>, using the full dataset in the time do¬ 
main (not using Taylor’s hypothesis). Fluctuations 
are defined as the deviation from the linear trend 
computed during a single tidal phase (6 hour inter¬ 
val). We compute the skewness and flatness of the 
fluctuations for each thermistor in order to assess 
the dependence on the height above the bottom. 
Results from upslope and downslope tidal phases 
are averaged separately, after normalisation by the 
standard deviation of the single tidal phase. 

Using the full dataset, skewness and flatness are 
determined by a combination of waves and turbu¬ 
lence that characterises a strongly stratified envi¬ 
ronment. In particular, there is a strong influ¬ 
ence from the waves at periods shorter than the 
tidal one, not removed by the linear detrending 
performed before computing the moments of the 
distribution. In an attempt to separate turbulence 
from waves, the data can be filtered. However, 
there is no obvious scale separation between waves 
and turbulence; this distinction may in fact be ill- 
defined. The fluctuations seem to be a continuum 
of waves of different degree of nonlinearity at fre¬ 
quencies lower than the (local) buoyancy frequency 
and turbulent motions at frequencies both below 
and, mainly, above N (see e.g. |Billant and Chomaz 
2001). Keeping in mind that the filtering can¬ 


as a function of the separation distance r (assum¬ 
ing homogeneity), usually along the horizontal di¬ 
rection (i.e. using spatial data exploiting Taylor’s 
hypothesis). Temperature increments are particu¬ 
larly useful when dealing with non-stationary time 
series (Davis et al. 1994). Increments enable one 


to naturally study turbulence and waves together 
without the need for filtering. 

Further insight into the turbulence dynamics 
can be gained by studying the behaviour of gen¬ 
eralised structure functions of temperature incre¬ 
ments (Warhaft |2000| ) 


Iq = lq{r) = (|Ar6»n . 


(3) 


This is done in section 13^ Generalised structure 
functions are a standard tool to characterise the 
intermittency of the velocity and scalar fields in 


a turbulent flow ( 

Frisch |1996 Shraiman and Sig-| 

|gia| 2000 

Warhaft 

2000). Here, we obviously fo- 


cus on the intermittency of the scalar (tempera¬ 
ture). According to the non-intermittent theory 
of Kolmogorov-Obukhov-Corrsin, the temperature 
field is expected to be fully self-similar within the 
inertial range, resulting in a structure function 


having a scaling behaviour of the form (Warhaft 

20001 ) 

with C,{q) = q/‘i^ 




(4) 


for every order q (see also the discussion of this 
in jCosta Frola et al.l |2014|). Note that a self¬ 


not separate, or define, these components, it still 
proves to be an instructive exercise, discussed in 
section [3^ We use a Butterworth filter removing 
frequencies below 1/5600 s“^, a frequency approx¬ 
imately equal to the lowest values of the buoyancy 


similar behaviour is compatible with a Gaussian 
distribution, but does not imply it. For g = 2, 
the structure function embeds the same informa¬ 
tion as the spectrum, with the slope C(2) = 2/3 of 
the structure function being the Fourier transform 
of the —5/3 slope of the fluctuations in wavenum¬ 
ber space. While for g = 2 the prediction of the 
theory is usually close to the values measured in 
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the laboratory, for g' > 2 there is an increasing dis¬ 
crepancy between the predicted g/3 and the ob¬ 
served C(^): fhe latter being systematically below 
the theoretical value, possibly reaching an asymp¬ 
totic value Coo: he. saturating, for sufficiently high 
orders {q ~ 10). This saturation is related to 
the presence of ramp-cliff structures in the scalar 
field, in connection with an accumulation of sharp 
gradients at small scales. A value of ({q) lower 
than expected is the consequence of an excess of 
these large gradients at small scales (the cliff struc¬ 
tures) and of a lack of large increments at large 
scales (the ramp structures). Evidence for satu¬ 
ration of moments has been reported in both pas¬ 
sive scalars (shear-dominated turbulence) and ac¬ 
tive scalars (convectively driven turbulence). How¬ 
ever, the estimated asymptotic value of ("oo is dif¬ 
ferent for the two cases. In particular, values of 
Coo ~ 1-4 are reported for passive scalars in grid 
turbulence ( Warhaft[ |2QQQ[ Celani et al. 2001), 
while smaller values (Ccxd ~ 0-8) are observed in con¬ 
vectively driven turbulence ( |Zhou and Xia 2002 


Celani et ah 


2002 ). 


In section inr we will consider “plus” and “mi¬ 
nus” increments, introduced by |Zhon and Xia 


( 2002| ) in the context of convectively driven scalar 
turbulence, and used here to better investigate the 
passive/active nature of temperature. They are 
defined as = (|A^6>| ± A^O) /2, and simply 

separate positive from negative increments in time 
(subscript r), measured at a fixed position. A plus 
(minus) increment is zero for all negative (posi¬ 
tive) increments, otherwise it is equal to the ab¬ 
solute value of the increment. A hot plume rising 
through a thermistor produces a steep increase in 
temperature with its leading cap, but a relatively 
gentle temperature decrease with its tail. This 
asymmetry between positive and negative tempera¬ 
ture increments should show up in the skewness of 
the plus and minus increment distributions, with 
higher skewness for the plus increments (note that 
plus and minus increments are positive definite). 
The laboratory results of Zhou and Xia ( |2QQ2 ) in¬ 
deed show this asymmetry in the convective region 
of a Rayleigh-Benard cell, for increments at in¬ 
ertial and dissipation separation distances. More 
generally, here we use plus and minus increments 
to study if ramp-cliff structures have a preferen¬ 
tial orientation. In our case, convective plumes, if 
present, are advected through the thermistor ar¬ 


ray according to Taylor’s hypothesis, rather than 
vertically rising through the array. Consequently, 
we have to consider spatial increments of temper¬ 
ature rather than time increments, and we do so 
along both the horizontal and along the vertical 
directions. Plus and minus increments are thus de¬ 
fined here as ArO^ = (|A^^| A ArO) /2 for incre¬ 
ments in the horizontal direction, and similarly as 
A;s6>^ = {\Az0\ ± AzO) /2 along the vertical. Here, 
AzO is defined by equationwith r substituted by 
the upward vertical coordinate. Given the rela¬ 
tively small number of sensors within each segment 
and the relatively low vertical resolution, AzO"^ is 
computed only at the thermistor spacing, i.e. 0.7 m. 


3 Results 


Velocities are higher during the upslope phase of 
the tide (see figureright panel, black triangles). 
Deviations of velocity from the vertical average ex¬ 
plain approximately 25% of the variance. Shear 
is computed at the 5 m vertical resolution of the 
ADCP, as the average of the magnitude of the ver¬ 
tical gradient of the horizontal (vector) velocity. 
Shear is similar during the two tidal phases, ap¬ 
proximately equal to 3.2 x 10“^ s“^ (as a mean of 
the instantaneous shear). Mean velocity is virtu¬ 
ally constant in the upper half of the mooring, and 
decreases by approximately 20% in the lower half. 
The temporal and spatial resolution of the ADCP 
is such that turbulent motions cannot be studied in 
any detail. Velocity has a strong component along 
the local isobath during both tidal phases. 

Figure hints at some important differences be¬ 
tween the two tidal phases. During the downs- 
lope phase, the sharp interface above the cold near¬ 
bottom layer is lower and thinner. Overturns at 
the sharp interface are also smaller than those ob¬ 
served during the upslope phase. Both these ele¬ 
ments suggest that turbulence is weaker during the 
downslope phase of the tide, as already observed by 
van Haren and Gostiaux (2012) and also by 


Mourn 


et al. (2004) in a shallower environment. This fact 


is confirmed by calculations on this dataset using 


similar methods to van Haren and Gostiaux (2012) 


(not shown). A peculiar feature of both panels in 
figure is that some of the turbulent structures, in 
particular in the upper half of the mooring, remind 
one of convective plumes. Particularly striking for 
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its size is the mushroom-like plume duriug the up- 
slope phase (black arrow iu the figure). 

3.1 Temperature horizontal 

wavenumber spectra 

As a startiug poiut iu the statistical characterisa- 
tiou of the temperature field, horizoutal waveuum- 
ber spectra of temperature fluctuatious compeu- 
sated by are showu iu figure The figure 

shows spectra averaged duriug the two tidal phases 
(upslope, dowuslope) aud iu each mooriug segmeut. 

Duriug the upslope tide (upper four spectra iu 
figure]^, au iuertial rauge is clearly visible iu all 
mooriug segmeuts as a horizoutal slope, approx¬ 
imately iu the k~^ rauge from 5 to 100 m. For 
segmeuts A aud B, there is also a hiut of a lo¬ 
cal miuimum at k~^ ~ 100 m, possibly the sigu 
of a separatiou betweeu the gravity-wave spectrum 
(lower waveuumbers, higher slope) aud the turbu- 
leuce baud. High-waveuumber roll-off of the spec¬ 
trum is visible above k~^ ~ 5 m, but the roll-off 
should be iuterpreted with care, siuce the 1 Hz sam- 
pliug rate of the thermistors is iusufhcieut to re¬ 
solve the dissipatiou scales (of both momeutum aud 
temperature) iu the most turbuleut portious of the 
dataset. Spectra from siugle tidal phases (6 hour 
slices of data), show that the high-waveuumber roll¬ 
off is uot always detected. Iu mauy cases, the spec¬ 
trum is flat up to almost the highest waveuumber, 
where the spectrum ofteu beuds upwards. This up¬ 
ward beud of the spectrum is due to the sigual ap- 
proachiug the uoise level of the measuriug device iu 
weaker stratificatiou aud less turbuleut regious. 

Iu the upper half of the mooriug, the spectra 
do uot reach the theoretical —5/3 slope, i.e. they 
slope dowuwards to the right iu figure We sug¬ 
gest that this is the result of averagiug together 
data regious haviug iuertial rauges spauuiug differ- 
eut scales. This effect is particularly stroug duriug 
the dowuslope phase, beiug less turbuleut thau the 
upslope oue. We couclude that the thermistors re¬ 
solve at least the lower-waveuumber part of the tur- 
buleuce iuertial rauge. Turbuleuce is strougly iuter- 
mitteut, iu particular duriug the dowuslope phase 
of the tide, as judged from the fact that the spectral 
slope is steeper thau —5/3. 

Similar couclusious could be drawu from the fre- 
queucy spectra (avoidiug the use of Taylor’s hy¬ 
pothesis). However, the spatial data have spectral 


slopes closer to the classical value of —5/3 duriug 
the upslope phase, aud also better highlight the 
differeuces betweeu the two tidal phases. 


3.2 Skewness and flatness of fluctu¬ 
ations 

Figure shows the skewuess aud flatuess of tem¬ 
perature fluctuatious, 6>, from the full dataset iu 
the time domaiu, as a fuuctiou of the height above 
the bottom aud of the tidal phase. The skewuess 
is always uegative aud small (figure left pauel). 
The most uegative skewuess appears iu the regiou 
above the uear-bottom cold layer, most likely iu 
couuectiou with mixiug of the cold layer itself; it de¬ 
cays towards zero further above. The flatuess (right 
pauel iu figure]^ is close to the Gaussiau value of 3, 
iu particular iu the upper half of the mooriug aud 
duriug the upslope phase of the tide. The Gaussiau 
value of the flatuess is clearly exceeded ouly close 
to the bottom duriug the dowuslope phase. 

As discussed iu sectiou |2.3[ fllteriug the data is 
a uaive attempt to separate turbuleuce from wave 
motious. Skewuess of the high-pass-filtered data 
(figure left pauel) is positive close to the bot¬ 
tom, while it is uegative further above. Deviatious 
from zero are larger thau iu the uufiltered data, as 
expected for a turbuleut sigual where uouliuear dy- 
uamics domiuate. Moviug towards the upper eud of 
the mooriug, skewuess slowly approaches zero. Af¬ 
ter fllteriug, skewuess is similar duriug the two tidal 
phases, iudicatiug that the differeuces observed iu 
figure 1^ cau mostly be ascribed to iuterual wave 
motious. The chauge of sigu of skewuess above 
the bottom is the siguature of the overturus that 
mix the iuterface above the cold uear-bottom layer. 
Skewuess has slightly larger uegative values duriug 
the upslope phase. 

The flatuess of the filtered dataset (figureright 
pauel) is much larger thau iu the full data; it is 
also larger duriug the dowuslope phase. Flatuess 
iucreases towards the bottom, iu particular iu the 
lowest 10 m of the mooriug. A flatuess larger that 
the Gaussiau value of 3, with expoueutial tails of 
the pdf as the oues fouud iu our data (uot showu), 
is a commou feature of turbuleut fluctuatious iu the 
iuertial aud dissipatiou rauges of both active scalars 


(Ghiug 1991) aud passive scalars iu the preseuce 


War haft 2000). 


of a meau gradieut (Shraimau aud Siggia 2000 
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From the comparison of filtered and unfiltered 
data, we conclude that even if the internal waves 
dominating the total signal seem highly nonlinear 
and continuously breaking, as suggested by figure 
their statistics are not too far from Gaussian. In the 
filtered data, on the other hand, turbulent fluctua¬ 
tions are strongly non-Gaussian. Furthermore, the 
skewness and flatness of the wave motions are af¬ 
fected by the phase of the tide more than the skew¬ 
ness and flatness of the turbulent motions. Both 
filtered and unfiltered data depend on the height 
above the bottom, in particular in segment A. 


3.3 Skewness and flatness of incre¬ 
ments 

We now consider the skewness and flatness of tem¬ 
perature increments, defined by equations and 
plotted in figure as a function of the separation 
distance r. Here, the values are presented as av¬ 
erages in the four segments. Note that A^O has a 
mean different from zero, since the tidal wave pro¬ 
duces a trend in 0; the mean increment is greater 
or smaller than zero depending on the tidal phase 
(downslope or upslope respectively). 

Skewness has a rather complex dependence on r. 
During both tidal phases, skewness increases from 
the bottom (A) to the top (D) of the mooring. For r 
smaller than 5 —10 m, skewness is higher during the 
upslope phase compared to the downslope phase, 
being positive only for segments G and D during 
the upslope phase. For separation distances larger 
than 10 m, skewness changes in opposite directions 
for the two tidal phases, towards more positive val¬ 
ues during the downslope phase and towards more 
negative values during the upslope phase. This be¬ 
haviour for larger scales is due to the opposite mean 
trend of temperature due to the tidal wave itself. 

In the inertial range, a skewness of the tempera¬ 
ture increments different from zero is considered the 


signature of ramp-cliff structures (Warhaft 2000). 
These ramp-cliff structures are most visible if in¬ 
crements are computed parallel to the mean scalar 
gradient (vertical), but they have a projection also 
on the horizontal plane. Here we focus on horizon¬ 
tal increments, owing to the limitations of the verti¬ 
cal resolution. The observed behaviour of skewness 
shows that the anisotropy of turbulence dynamics 
(as measured by skewness) depends on the phase 
of the tide and on the height above the bottom at 


the smallest scales resolved in the dataset. If skew¬ 
ness is determined by small-scale convective plumes 
(producing structures that can be assimilated to 
ramp-cliff ones, Zhou and Xia, 2002), the observed 
differences in skewness could be caused by the dom¬ 
inance of either warm or cold convective plumes, 
depending on the tidal phases. Thorpe et al. (1991) 
suggested that negative skewness is associated with 
convection, using observations from near-surface 
and near-bottom oceanic boundary layers. Here, 
the clearest signs of convection are found in seg¬ 
ments G and D during the upslope tide (as will be 
discussed in sections 3.4 and 3.5), having positive 
skewness at small separations. The difference be¬ 
tween our results and those of Thorpe et al. (1991) 
could be caused by the prevalence here of warm 
rather than cold plumes, a hypothesis that is dis¬ 
cussed in more detail below. More generally, we 
suggest that the skewness of the temperature in¬ 
crements is not a good indicator of the presence or 
absence of convective activity. 

Flatness decreases for increasing separation dis¬ 
tances during both tidal phases, approaching with¬ 
out actually reaching the Gaussian value of 3 at the 
largest scales considered. A broad range of power- 
law decay is evident, as often observed within the 
inertial range in the laboratory (e.g. figure 3 in 
Warhaft 2000). Flatness is expected to saturate 


for r comparable to the dissipation scale, but we 
barely see a hint of such a saturation at the small¬ 
est r, confirming that the dissipation scale is not 
resolved here. Flatness is very similar throughout 
the mooring during the downslope phase, while it 
is higher away from the bottom during the upslope 
phase. 

In general, skewness and flatness strongly depend 
on r. Skewness in particular has a rich behaviour 
within the turbulence inertial range identified in the 
spectra in figure However, variations in skew¬ 
ness are generally smooth, and overall these quan¬ 
tities do not provide any clear suggestion of a “scale 
break” between turbulence and wave dynamics. In 
the next section, we will show that, on the other 
hand, higher-order moments give a clear indication 
of such a scale break. 


3.4 High-order moments 

We extend the description of the statistics of the 
temperature increments to higher-order moments. 
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using in particular the generalised structure func¬ 
tion defined in [3l In order to measure the con¬ 
vergence of the computation of these higher-order 
moments, we consider the quantity 

max|Ar6»|® 

M%ir) ’ ^ ^ 

where M is the number of data points used to 
compute the average 7 ^. This ratio is the relative 
contribution of the largest resolved increment (at 
scale r) to the generalised structure function. A 
well-converged generalised structure function im¬ 
plies a small value of this ratio, since the incre¬ 
ments most relevant for the order q are then well 
resolved, i.e. they are smaller than max|A^^|. In 
order to improve the convergence of 7 ^ for higher q, 
we consider the statistics averaged in two vertical 
subsets of the mooring instead of four as done in the 
previous sections. We consider the lower half (seg¬ 
ments A and B combined together) and the upper 
half (segments C and D combined together). 

Based on the ratio jq is well converged at 
scales larger than 1 m at all orders (ratio 0 ( 10 “^)), 
with the partial exception of the highest orders in 
the lower half of the mooring during the upslope 
phase (not shown). For scales smaller than 1 m, 
the convergence is worse in particular during the 
upslope phase, and the orders higher than q ^ 6 
are unreliable. For this reason, we will focus on 
scales larger than 1 m. 

Figure shows, in log-log space, the values of 
7 g as a function of r, for q ranging from 1 to 10 , 
in the lower and upper halves of the mooring and 
during the two tidal phases separately. The scal¬ 
ing exponents Q (equation are estimated in the 
turbulence scaling range and are shown in figure 
as a function of the order q. Considering the lower 
half (A-l-B panels in figure]^, 72 has a scaling ex¬ 
ponent close to 2/3 over almost two decades during 
the upslope phase. A scaling range with exponent 
1.4 (dotted lines) up to r ^ 100 m is approached 
for high q in the case reported in the lower left 
panel of figure Figure however, shows that 
the scaling exponents do not saturate at this value, 
but keep increasing with q. During the downslope 
phase (figure lower right panel), 72 has a scal¬ 
ing exponent higher than 2/3, and 7 ^ has an ex¬ 
ponent larger than 1.4 at the highest orders com¬ 
puted. A clear break of the scaling behaviour is 
visible at r ~ 100 m for g' > 5, in particular during 


the downslope phase. It is unlikely that this is a 
consequence of the violation of Taylor’s hypothesis, 
because the scale of the break is smaller during the 
upslope phase, having larger mean velocity. We 
would expect on the other hand the scale of the 
break to be proportional to the mean velocity, if it 
were caused by a breakdown of Taylor’s hypothesis. 
Higher-order moments thus identify two separate 
ranges, which could be interpreted as turbulence- 
dominated for r smaller than approximately 100 m 
and wave-dominated for larger scales. 

Considering now the upper half of the mooring 
(C+D panels in figure [^, 72 has a scaling exponent 
greater than 2/3, markedly so during the downslope 
phase. The scaling behaviour at high q extends 
up to scales not as large as in segments A+B (30- 
50 m). Wave motions may thus reach smaller scales 
in the upper half of the mooring, and turbulent mo¬ 
tions may be limited to a narrower range. The value 
of Cg for the highest g’s is well below 1.4 during the 
upslope phase (see figure |^, while it is higher than 
1.4 during the downslope one. Considering that 
72 has a slope in excess of 2/3, we hypothesise that 
also the slopes for higher q are similarly biased with 
respect to controlled environment results. If this is 
actually the case, the asymptotic exponent during 
the upslope phase would in fact approach the value 
of 0.8 measured for an active scalar in convective 


turbulence (Zhou and Xia 2002 ). 


Considering the scaling exponents in figure 
more generally, we see strong deviations from the 
non-intermittent = g/3 prediction, marked by 

the dashed line in the figure. The strongest devia¬ 
tions are observed during the upslope tidal phase, 
in particular in the upper half of the mooring 
(C-j-D). Saturation of the moments is clearly vis¬ 
ible in this part of the dataset, already for g = 5. 
These higher-order moments thus show that the 
tidal phase actually has a strong impact on inter- 
mittency, and that the turbulence behaviour is far 
from universal in the dataset. 


3.5 Plus and minus increments 

In the previous section, we suggested that the be¬ 
haviour of temperature is similar to that of a pas¬ 
sive scalar in the lower half of the mooring, while in 
the higher part the statistics deviate towards those 
of active scalars, in particular during the upslope 
phase. To further investigate this issue, we con- 
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sider the skewness of plus and minus increments, 
introduced in section [231 The skewness of the dis¬ 
tributions of and are shown in figure 

as a function of the separation distance and for 
the four segments of the mooring separately. 

During the downslope phase, the skewness of the 
vertical minus increments, fis [AzO~) (large, filled 
red triangles), is greater than //s (large, 

filled red squares) in all segments. During the ups- 
lope phase of the tide this is also true, but the dif¬ 
ference is much smaller in the lower segments that 
in upper ones. Considering the horizontal incre¬ 
ments during the downslope phase of the tide, the 
skewness of the minus increments (small, red filled 
triangles) is larger at all scales and throughout the 
mooring. In the upslope phase closer to the bot¬ 
tom (segments A and B), the skewness of plus and 
minus horizontal increments are the same, or close, 
with the exception of larger scales. On the other 
hand, plus skewness (small, blue open squares) is 
larger in segments C and D at all scales. 


The difference in skewness between plus and mi¬ 
nus increments is of order 1, consistent with that 
reported in the laboratory by | Zhou and Xia (2002). 
On the other hand, laboratory data show that 
fis (At-6 >+) and /is (At-6 >“) converge to the same 
value for large r. This does not happen here, where 
an increase of the difference is observed for r larger 
than 10-100 m. We should not conclude from this 
that plumes of all sizes are present, but rather that 
wave motions are also associated with asymmetric 
temperature gradients. We also cannot rule out 
that contamination between the two tidal phases 
and the use of the Taylor’s hypothesis may play a 
role at the largest scales. 


These results are discussed with the help of the 
sketch in figure We consider two “plumes”, 
which we regard here generically as an asymmetric 
temperature anomaly, with a sharp leading front 
and a smoother tail. In the upper part of the 
figure, the isotherms of two of these plumes are 
sketched in the spatial domain: (1) on the left a 
hot plume with an upper sharper front; and (2) on 
the right, a cold plume with a lower sharper front. 
Two transects are taken through each of the two 
plumes, one vertical (green dashed arrow) and one 
horizontal (red continuous arrow). The resulting 
temperature transects are sketched in the two pan¬ 
els below, with green dashed and red continuous 


lines for vertical and horizontal directions, respec¬ 
tively. Consider the hot plume on the left. Its 
leading front is sharper than its tail; the vertical 
transect thus shows a slower increase of the tem¬ 
perature and a faster decrease to the mean value. 
This structure of the plume shows up as larger val¬ 
ues of skewness for minus rather than plus verti¬ 
cal increments, up to separations comparable with 
the size of the plume. The same reasoning applies 
to the horizontal direction and for a plume with 
an inclination to the mean flow given in the fig¬ 
ure. The result is a larger skewness for plus than 
for minus horizontal increments. Note that a cold 
plume with opposite tilt would produce different re¬ 
sults for the vertical increments, so the asymmetry 
is characterised unambiguously. This is consistent 
with the estimates in segments C and D during 
the upslope tidal phase. A similar reasoning can 
be applied to the cold plume in figure suggest¬ 
ing that in this case /is (^Az0~^) < /is (^Az0~) and 
^3 {Ar0+) <113 (Ard-). 

We can now interpret the results from figure as 
follows. During the upslope phase, warm plume- 
like temperature anomalies are dominant in seg¬ 
ments C and D, well above the bottom. Dur¬ 
ing the downslope phase, cold plume-like anoma¬ 
lies are instead dominant throughout the mooring. 
The asymmetry of temperature anomalies has both 
a vertical and a horizontal component. In other 
words, plumes also have a preferential horizontal 
orientation, pointing towards the direction of the 
mean flow, as shown in figure On the other 
hand, in segments A and B during the upslope tide, 
temperature anomalies are more symmetric, in par¬ 
ticular along the horizontal direction. Overall, we 
do not have any clear indication of the size of the 
plumes. We can only speculate that their typical 
size is at least of 0(1 m), i.e. at least comparable 
to the values of r for which skewness departs from 
its saturated, small-scale value. The large plume 
of figure is thus an extreme case, and scanning 
through the temperature record confirms this. 


4 Discussion and conclusions 

We described the temperature statistics using ob¬ 
servations from a moored thermistor array, de¬ 
ployed in the deep ocean above a sloping bot¬ 
tom. The dynamics at the mooring location are 
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strongly affected by the semidiurnal tidal wave, 
and we thus considered separately the statistics for 
two tidal phases, cooling (upslope) and warming 
(downslope). The thermistor array does not reach 
the bottom, with the lowest sensor approximately 
5 m above the bottom, and a bottom mixed layer, 
if present, is not visible in the data. The results 
are one of the few studies on the high-order scalar 
statistics in stably, strongly stratified fluids. 

In the dataset, temperature fluctuations and 
temperature increments are strongly non-Gaussian 
at (horizontal) scales smaller than 10-100 m. At 
larger scales, deviations from Gaussian are smaller, 
but there generally is still a scaling behaviour, with 
a higher value of the scaling exponent Q. The sur¬ 
prising result here is that high-order moments show 
an abrupt transition between the two regimes (see 
in particular panels A-^B up and down, G+D down 
in figure [^. The length scale of the transition 
is of order 100 m, but the actual value is state- 
dependent, with significant variations depending 
on the tidal phase and on the height above the 
bottom. The 100 m scale is equal to the buoy¬ 
ancy length obtained as the ratio of the velocity 
scale (10“^ ms“^; see figureand the buoyancy 
frequency (10“^ s“^, see table |^, and is approx¬ 
imately 100 times larger than the Ozmidov scale 
(estimated from the Thorpe scale, averaged over 
the whole dataset, computed in |Gimatoribns et al.[ 
2014| ). Since the layer thickness scales with the 
buoyancy length in a stratified turbulent environ¬ 
ment (Billant and Ghomaz 2001), we hypothesise 


(state-dependent) scale separation between waves 
and turbulence may have important practical appli¬ 


cations, e.g. for large eddy simulations (Gosta Frola 


et al. 2014). It is also noteworthy that, assum¬ 


ing we can indeed identify large-scale motions with 
waves, deviations from Gaussian are small even if 
there is ample evidence of nonlinearity and break¬ 
ing. 

We also considered generalised structure func¬ 
tions and skewness of plus/minus increments with 
respect to the passive/active nature of temperature 
dynamics. There is a clear indication of a passive 
scalar behaviour during the upslope phase close to 
the bottom; shear is strong enough to dominate 
the buoyancy forces. On the other hand, in the up¬ 
per half of the mooring, temperature statistics ap¬ 
proach those of an active scalar during the upslope 
phase. The difference in skewness of plus/minus 
increments is consistent with the presence of warm 
convective plumes. The numerical simulations of 


Slinn and Levine (2003) and Gayen and Sarkar 
]2011 a|b ) suggested that tidal motions above slopes 


that motions at larger horizontal scales are less non¬ 
linear, i.e. more wave-like, as overturns are lim¬ 
ited by the layer thickness. In this view, we ar¬ 
gue that this break in the scaling behaviour could 
be a characteristic feature of stratified turbulence. 
The fact that a sharp break is visible only in the 
high-order moments, at least in this dataset, could 
result from the fact that the buoyancy length is 
relevant only for turbulent motions. On the other 
hand, internal waves are scale-free, at least ac¬ 
cording to the dispersion relation from linear the¬ 
ory, and thus have scales both smaller and larger 
than the break. Since waves are less effective than 
turbulent motions at producing sharp gradients, 
their signature on higher-order moments is smaller 
than for turbulent fluctuations, and thus smooth 
the transition between the two regimes only in the 
lower-order moments. The possibility to define a 


can produce strong convection and mixing, in par¬ 
ticular during the upslope tidal phase. Our results 
provide support for their conclusions. 

During the downslope phase, results are more 
ambiguous. The difference in skewness of 
plus/minus increments indicates an asymmetry of 
temperature anomalies consistent with the presence 
of cold plumes throughout the mooring. However, 
the scaling exponents of the generalised structure 
function exceed those of passive scalars (figure]^. 
The fact that in this case low-order moments exceed 
the q/3 prediction suggests, however, that this less 
turbulent tidal phase is strongly inhomogeneous, 
and a comparison with laboratory results is more 
difficult. We hypothesise that the observed scaling 
behaviour is typical of strongly stratified, sporadi¬ 
cally turbulent flows. In this respect, it would be 
useful to compare our results with numerical sim¬ 
ulations similar to those reported in |Rorai et al. 


(2014), if scalar statistics can be computed from 


them. 

In conclusion, we want to return to the initial, 
broader, issue of mixing in geophysical flows. 
Overall, the indication from our results is that 
temperature statistics are highly dependent on the 
state of the flow, even if stratification is always 
present. Differences in the statistics, in particular 
in relation to convective motions, are very likely 
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to be connected to changes in mixing properties, 
in particular mixing efficiency. A clearer under¬ 
standing of mixing in a stably stratified flow is 
an important open question with implications 
in the fields of oceanography and climate. We 
hope that the present results can help in this effort. 

The authors would like to thank the crew of 
RV Pelagia for enabling the deployment and re¬ 
covery of the instruments, and all the technicians, 
in particular Martin Laan, indispensable for the 
success of the measurements. The authors would 
like to thank the referees for their useful and de¬ 
tailed comments on the previous versions of the 
manuscript. 
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Latitude 

Longitude 

Deepest thermistor depth 
Deepest thermistor h.a.b. 
Bottom slope 
Number of thermistors 
Thermistor vertical spacing 
Total length of array 
Total length of cable 
Deployment 
Recovery 


36° 58.885'N 
13° 45.523'W 
2205 m 
5 m 
9.4° 

144 
0.7 m 
100.1 m 
205 m 

13 Apr 2013 
12 Aug 2013 


Table 1: Details of the mooring from which the data is obtained, h.a.b. is the height above the bottom. 


Segment name 

A 

B 

c 

D 

h.a.b. range 
Thermistor number 

5.0-30.0 m 
1-36 

30.0-55.0 m 
37-72 

55.0-80.1 m 

73-108 

80.1-105.1 m 
109-144 


Table 2: Definition of the segments used in the analysis of the data, h.a.b. is the height above the 
bottom. 


Quantity 

Tidal phase 

Segments 

N 

A 

B 

c 

D 

all 

Mean 

Both 

1.46 

1.23 

1.26 

1.33 

1.32 


Up 

1.43 

1.28 

1.29 

1.33 

1.33 


Down 

1.49 

1.17 

1.22 

1.34 

1.30 

Percentile 1 

Both 

0.21 

0.18 

0.19 

0.21 

0.19 


Up 

0.21 

0.19 

0.19 

0.21 

0.20 


Down 

0.20 

0.17 

0.18 

0.21 

0.19 

Percentile 99 

Both 

6.15 

5.12 

5.06 

5.10 

5.40 


Up 

5.98 

5.42 

5.32 

5.17 

5.48 


Down 

6.29 

4.79 

4.79 

5.03 

5.31 


Table 3: Mean, 1th and 99th percentiles of buoyancy frequency (from reordered 25 s average profiles of 
temperature) in the four different segments of the mooring and for the full dataset. The results are in 
units of 10“^ s“^. 
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/ ^ [hour] V, (v ), f'[m s 


Figure 1: Spectra (left) and pdf’s (right) of total velocity v (black), mean velocity {v) (red, see text) and 
velocity fluctuations v' (light blue). In the left panel, the thick dashed line shows the M 2 semidiurnal 
frequency. On the right panel, the lines show the pdf for all the data, while triangles (circles) show the 
pdf computed for the upslope (downslope) tidal phase alone. Only velocities measured within the depth 
range of the thermistor array are used in the figure. 



Figure 2: Examples from the dataset during a downslope (down) and upslope (up) tidal phase. The 
dashed lines mark the limits of the four segments of the mooring (A, B, C, D) considered separately 
throughout the paper. Colour scale is different in the two panels. The lowest thermistor is located 
approximately 5 m above the bottom. The vertical scale is in meters above the bottom. The black 
arrow points to the large plume mentioned in the text. 
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Figure 3: Average wavenumber spectra of temperature fluctuations, during the two tidal phases and in 
the four segments of the mooring, compensated by An arbitrary vertical shift is applied to the 

spectra to better distinguish them. 
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Figure 4: Skewness (left) and flatness (right) of temperature fluctuations, plotted as a function of depth 
for every other sensor in the string, and computed separately for different tidal phases. Results from the 
unfiltered time series (without using Taylor’s hypothesis). The thick black lines mark the values of a 
Gaussian distribution. Dashed lines and letters refer to the different segments of the mooring (see table 
[^. Open blue triangles refer to the upslope phase, filled red circles refer to the downslope phase. 
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Figure 5: As figure]^ but for filtered data (see text). As a consequence of the filtering, the estimate of 
the flatness is noisier than in figure 
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Figure 6: Skewness (/is (A^^)) and flatness (/i 4 (A^^)) of temperature increments for the four segments 
of the mooring, computed for upslope and downslope phase as a function of the separation r. The 
black lines mark the values for a Gaussian distribution. Note that the vertical axis of the right panel 
is logarithmic. Different colours refer to different tidal phases (open blue for upslope, filled red for 
downslope), different symbols to different segments of the mooring (see legend). 
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Figure 7: Generalised structure function of temperature increments 7 ^. In the upper (lower) panels 
the results for the segments C, D (A, B) combined together are shown. The left panels refer to the 
upslope phase of the tide, while the right panels refer to the downslope phase. Different orders q are 
plotted, increasing from top to bottom in each panel, with colours changing with q (see legend). The 
black dashed line near the top is a guide for the eye with a slope of 2/3, the theoretical slope for the 
structure function of order 2. The dotted black lines in the lower part of each panel have a slope of 1.4, 


the asymptotic value of the slope in the inertial range for a passive scalar in grid turbulence (Warhaft 


2000 ). 
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Figure 8: Exponents C,q of the generalised structure functions shown in figure as a function of the 
order g, estimated by a least square fit in the scaling range of turbulence. The dashed line is the 
C <7 = slope for the non-intermittent case. The dotted horizontal line marks C = 1.4, the asymptotic 
value approached in grid turbulence for large q. Filled red (open blue) symbols refer to the downslope 
(upslope) tidal phase. Circles (stars) refer to the lower (upper) half of the mooring. 
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Figure 9: Skewness of plus and minus increments as a function of the separation distance r. The four 
panels refer to the A-D segments of the mooring. Colours refer to different tidal phases (open blue 
for upslope, filled red for downslope), different symbols to plus and minus increments (squares for plus, 
triangles for minus). The skewness of both horizontal increments //s (small symbols) and vertical 

increments /is (A;^^^) (large symbols) are reported (see legend). 
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Figure 10: 


Sketch of the signature of convective plumes found in the dataset, see text for an explanation. 
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